Adiabatic Creation of Atomic Squeezing in Dark States vs. Decoherences 
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We study the multipartite correlations of the multi-atom dark states, which are characterized by 
the atomic squeezing beyond the pairwise entanglement. It is shown that, in the photon storage 
process with atomic ensemble via electromagnetically induced transparency (EIT) mechanism, the 
atomic squeezing and the pairwise entanglement can be created by adiabatically manipulating the 
Rabi frequency of the classical light field on the atomic ensemble. We also consider the sudden 
death for the atomic squeezing and the pairwise entanglement under various decoherence channels. 
An optimal time for generating the greatest atomic squeezing and pairwise entanglement is obtained 
by studying in details the competition between the adiabatic creation of quantum correlation in the 
atomic ensemble and the decoherence that we describe with three typical decoherence channels. 
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I. INTRODUCTION 

Atomic ensemble can serve as a quantum memory for 
storing the quantum information of photons where 
the memory elements can be described as the quasi-spin 
wave excitations in atomic ensemble [5], or Dicke type 
collective states [6]. Recent experiments have demon- 
strated that such storable and quantum memory could 
have long lifetime for the use in long-distance quantum 
communication [7\. When using the magnetically insen- 
sitive clock transition in atomic rubidium confined in a 
one-dimensional optical lattice, quantum memory life- 
time can exceed 3 ms [8]. 

The relevant quantum storage process could be im- 
plemented through the many-particle enhancement of 
the absorption cross section with the adiabatic-passage 
techniques [9.] . Through the electromagnetically induced 
transparency (EIT) mechanism, the adiabatic manipula- 
tion utilizes the dark polariton \dn{0)) (a photon-dressed 
collective atomic state). With adiabatically changing 
parameter 6 G [0,7r/2], the quantum state of photon 
\P) = X]^Cn|n) (a superposition of photon Fock states 
\n)) can be adiabatically converted to a state of the col- 
lective atom excitation \M) = Xln^^l^^)' which is a 
superposition of the collective atomic states \Mn) with 
the same coefficients as that in the photon state |P). 
Mathematically, this coherent conversion of the photon 
state is an associated mapping 



IP) ® \d') 



;c„|d„(0)) ^ ^c„|d„(|)) = \p') ® \M), 

(1) 

where \d^) and \p^) are the initial state of the memory 
and the final state of photon respectively. Obviously, 
here only used are the macroscopically coherent proper- 
ties of each collective excitation, which is described by 
individual collective state |M^), rather than the inter- 
nal entanglement and quantum correlation in the single 
collective state \Mn) as well as \dn{0)). 

In this paper we will pay attention to the latter by 
considering the multipartite correlation measured by spin 



squeezing in |Mn), which is obtained from \dn{0)) by adi- 
abatic manipulation for the dark state \dn{0)) as 



\n) |0) = Km 



|O)0|M,). (2) 



We study in details the dynamic competition between the 
adiabatic creation of quantum correlation in the atomic 
ensemble and the decoherence that we describe with three 
typical decoherence channels. There are two time scales 
ti and t2 depicting such competition, where ti repre- 
sents the adiabatic time limited by the adiabatic condi- 
tions and t2 represents the decoherence time. Actually, if 
we regard the atomic ensemble as a macroscopic object 
in large N limit, quantum decoherence scenario seems 
to result in the quantum disappearance due to noisy 
ambient environments [10] or the internal random mo- 
tions fill, [T2|. A surprising discovery about the pairwise 
entanglement is that it suddenly die due to influence of 
environment pisl, fl3|. In this paper, we also consider this 
sudden death phenomena of the pairwise entanglement 
as well as the atomic squeezing for quantum correlation 
for more than two particles. We compare the time scales 
of various sudden death processes with the speed of adi- 
abatic manipulations for creating such squeezing as well 
as the pairwise entanglement. 

Actually, creating spin squeezing in a single Dicke like 
state is the crucial step for the precision measurements 
based on many-atom spectroscopy. With such controlled 
production of complex entangled states of matter and 
light, squeezing the fluctuations via entanglement be- 
tween 2-level atoms can improve the precision of sens- 
ing, clocks, metrology, and spectroscopy [15]. This paper 
seems to provide a simple scheme for adiabatic creation 
of atomic squeezing. This adiabatic passage manipula- 
tion seems feasible in creating the spin squeezing, but it 
requires an ideal technique of single photon source, which 
can initially produce the single Fock states |n). Therefore, 
our present theoretical proposal may meet the difficulty 
for physical implementation based on the existing tech- 
nology. 

This paper is arranged as follows. In Sec. [Ill we study 
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FIG. 1: (Color online) Schematic illustration of the atomic 
ensemble interacting with two light fields: the classical light 
field (denoted by the blue arrow) with Rabi frequency Q and 
the wave vector i^me and the quantum light field (denoted 
by the red arrow) with the coupling constant g and the wave 
vector i^ge. Each A— type atom confined in the rectangular 
container has the same excited state |e), the relative ground 
state 1^) and the metastable state |m). 



the multipartite correlations characterized by the atomic 
squeezing beyond the pairwise entanglement in the dark 
states. In Sec. |TT1 we study the relationship between the 
concurrence and the spin squeezing under adiabatic ma- 
nipulation, and indicate that the photon statistics varies 
in the same way of the spin squeezing. In Sec. IIV| the 
sudden death of the atomic squeezing and the pairwise 
entanglement are considered by introducing three typical 
decoherence channels. The optimal times for generating 
the atomic squeezing and the pairwise entanglement are 
obtained in Sec. (Vl We conclude in Sec. (Vl 



II. ADIABATICALLY CREATING SQUEEZING 
IN DARK STATES 

A. Dark states 

As shown in Fig. [H the model we consider consists 
of atomic ensemble with TV A— type atoms, which have 
the excited state |e), the relative ground state \g) and 
the metastable state |m). For convenience we assume all 
the atoms have the same energy spacings for the above 
three atomic states. For the cold atomic ensemble, the 
Doppler width broadening of the atoms can be depressed 
effectively. These atoms interact with two single-mode 
light fields: one is quantized radiation mode (with cou- 
pling constant g and annihilation operator a) to lead the 



approximately resonant transition from |e) to l^f); the 
other is an exact classical field with Rabi frequency to 
lead the transition from |e) to |m). The quantum dy- 
namics of the total system is described by the following 
Hamiltonian in the interaction picture ^] 



N 

H = gaY^ exp (zKge • j) a^^^ 
j=i 

N 

n ^ exp (iKme • j) CTme + ^'C' 



(3) 



where Kge and K^e are, respectively, the wave vectors 
of the quantum and the classical light fields. We have in- 
troduced the quasi-spin operators cr^^ = (a, P = 
e, ^, m) for a ^ /3 in above Hamiltonian to describe the 
transitions among the levels of |e) , \g) and |m) . 

For EIT, to find the degenerate class of the eigenstates 
\dn (<9)) satisfying 

H\dn{0))=0, (4) 

we introduce the dark state polariton operator 

sin6>, (5) 

which obviously mix the electromagnetic field with cre- 
ation operator and the quasi-spin wave with corre- 
sponding operator 



1 



N 



,^^exp {iK-i)ai^g. 
j=i 



(6) 



The angle tan^ = g^/njO. can be adiabatically manip- 
ulated in the quantum storage process, and K = Kge — 
Kme is the wave vector difference between the quantum 
and the classical fields. We construct the degenerate class 
of the eigenstates \dn (0)) as 



(7) 



where the normalization constant 
A{n,N,0) -- 



N\ sin^^ 20 
[ iFi(-n;7V-n + l;-7Vcot2|9)]" 



(8) 



reduces to 1 when the atom number N is much larger 
than n. Here, n is the excitation number of the atomic 
ensemble, the vacuum state 



10) 



|0) ® \i) 



(9) 



is the direct product of the photonic vacuum state 
|0) and the quasi-spin ground state ||) = Ilj^i • 
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Here, iFi(a; 6; c) is the Kummer hypergeometric func- 
tion [1^. Eq. (jl]) indicates that \dn {0)) is totahy can- 
cehed by the interaction Hamiltonian, and thus is cahed 
a dark state or a dark state polariton. 

Obviously, the dark state \dn (6)) is photon Fock state 
\n) when ^ = and quasi-spin wave state when = 
7r/2. In the low excitation limit as n <C A/", CMs an 
approximately bosonic operator as well as (6). 

With the help of the dark states, by adiabatically 
manipulating the classical field the quantum informa- 
tion of the photon with quantum state |P) = ^^Cn \n) 
can be coherently converted into the atomic ensemble as 
\M) = Cn |M^), where \n) is the Fock state of photon 
and 



\Mn) 



{-If 



N^{N -n)\ 



N\n\ 



[Ct]"|0) (10) 



is the collective atomic state. This coherent conversion 
is an associated mapping 



|p) ® I ;) 



.Mn(0))^^C„|d„(|)) = |0)$5|M), 

(11) 

which only use the macroscopically coherent properties 
of each collective excitation described by individual col- 
lective state \Mn) . We notice that |M^) is actual multi- 
atom state with multipartite quantum correlations. For 
two-particle case, we call it internal entanglement. 

The internal entanglement and the quantum correla- 
tion in the single collective state \Mn) as well as \dniO)) 
also play important role in quantum storage process. 
Now we consider the quantum correlation measured by 
spin squeezing of those quasi-spins in IM^), which is ob- 
tained from \dn{0)) by adiabatic manipulation 



|n) ® |0) = |d„(0)) ^ M„(-)) 



|0) ® |M„). (12) 



B. Atomic squeezing in dark states without 
decoherence 



The above dark state mixes the electromagnetic field 
and the quasi-spins defined by the two levels \g) and |m) 
of the atomic ensemble. The collective operators 



1 ^ 



(13) 



can be generally used to demonstrate the exchange sym- 
metry and dynamic properties of the quasi-spin system. 
To calculate the averages of the collective operators, the 
reduced density matrix of the quasi-spins is obtained by 
tracing over the degree of freedom of the photons as 



pl = Tr^\dn{e)){dnm or 



E 

/e=0 



n! (cos^)^^ (sin6>) 



\2n-2k 



\A(n,N,e)\^ N\ 



k\{n-k)\ 



U 

where 



TV 

y 



^ 2 



N 

y 



N^{N -n^k)V 



(14) 



(15) 



is a symmetric Dicke state and the collective operator 
J+ = VnU^C^U (16) 

is a unitary transformation of the quai-spin wave opera- 
tor with the unitary matrix 



U = exp 



N 



3=1 



. H + 1) 



(17) 



Since only is the z— component quasi-spin operator 
contained in the unitary transformation, which will never 
change the quantum number of the symmetric Dicke 
bases, the density matrix in Eq. (p!4|) only has the diag- 
onal elements. Thus the averages of Jx and Jy vanish as 
{Jx) = {Jy) = 0- III this situation, only the 2:— component 
collective operator survives as 



(18) 



where {Jz)\q=t^i2 = ^ — N/2 is mean component spin 
of the Dicke state \n — N/2) and 



nN cot^ 6>iFi (1 - n; A/" - n + 2; -N cot^ 



{N -n^l)iFi{-n;N 



l;-A/'cot2 0) 



(19) 

is the deviation of the 2:— component spin. The fact that 
{Jz) does not vanish means that the mean spin is along 
the direction. Since SJz vanishes when = 7r/2, its 
value can measure the mixture of the electromagnetic 
field and the quasi-spin wave. 

The spin squeezing has several measures, and we only 
list three typical and related ones as follows: 



4(J2) 

S3 



(20a) 
(20b) 
(20c) 



In Eq. (j20ap . the minimization is over all the direction 
denoted by ±, which is perpendicular to the mean spin 
direction (J) / (J^) . For the quasi-spins in the dark state. 
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the and ^2 actually measure the spin squeezing of the 
atomic ensemble in the x — y plane. In Eq. 
the minimal eigenvalue of matrix 



r = (7V-l)7 + C, 



(21) 



where 



'kl 



{Jk){Ji) {k,le{x,y,z}) (22) 



is the covariance matrix and 



^kl = -^{JlJk + JkJl) 



(23) 



is the global correlation matrix. Parameters ^1,^25 
^3 were defined by Kitagawa and Ueda, Wineland et al., 
and Toth et al., respectively. If ^2 < ^ (^3 < ^)-> ^Pi^ 
squeezing occurs, and we can safely say that the multi- 
partite state is entangled [H, [l^l • Although we cannot 
say that the squeezed state according to < 1 is entan- 
gled, it is indeed closely related to quantum entanglement 

For the quasi- spins in the dark state, since the averages 
{JxJz) = {JyJz) = vanish as the same reason of the 
vanishing of {Jx) and (Jy), the and ^3 are simplified 

as M.m 



e. = iA{Ji^Ji)-\m\) 



N 

2 min{^?,c^2| 



2_ ' 

AT 



(24a) 
(24b) 




FIG. 2: (Color online) The squeezing parameter Q\ versus 
the excitation number n for ^ = (blue dashed line), — n /A 
(orange dotdashed line) and — n/2 (red solid line). The 
total number of the atoms N are chosen as 20. The horizontal 
purple dotted line as the base line represents the upper limit of 
the spin-squeezing. The dark state is always squeezed except 
without excitation (n = 0) or fully excited {n — N) . 



and the straightforward calculation gives 

5Jl = N cot^ e[2n + N cot^ 9 
{N +l){n + Ncoe 9) 



(iV - n + 1) 



T(n,N, 



(28a) 
(28b) 

(28c) 



where 



where 



A/-2 



N{^J,f + {J,f 



(25) 



characterizes the spin squeezing along the z— direction. 
Because {J-) = for the quasi-spins in the dark state, 
^2 and ^3 are always greater than 1 and thus there is no 



spin squeezing of the atomic ensemble in the x - 
In the following discussions, we only use 
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Ar2 



- y plane. 



(26) 



to characterize the 2:— component spin squeezing of the 
atomic ensemble in the dark states. 

We consider the special case that when the wave vector 
difference K between the quantum and the classical field 
is zero. Since /7|k=o identity matrix, the quai-spin 
wave operator in Eq. (p!6|) actually is proportional to the 
collective operator and the reduced density matrix in Eq. 
(p!4|) is spanned by the symmetric Dicke states with total 
quasi-spin J = N/2. Thus 



(27) 



N 
"2" 



(29) 



is mean z— component spin square of the Dicke state 
|n - N/2) and 

^("'^'^^= ,F,(-n;iV-n+l;-iVcot^^) - ^'^^ 

For convenience we define equivalent squeezing parame- 
ter El 



C^=max{l-CiO}, 



(31) 



which characterizes the spin squeezing when < C3 ^ 1- 



C. Numerical pictures of the atomic squeezing in 
dark states 

The squeezing parameters of the dark state \dn (0)) 
without wave vector difference (K = 0) are illustrated in 
Fig. [21 [3] and [H In all three figures, the total number of 
the atoms N are chosen as 20. Since the reduced density 
matrix in Eq. (p!4|) is symmetric when is replaced by 
6> + TT, the range of the 6 is chosen as [0, tt). Fig. [2] shows 



FIG. 3: (Color online) The polar plot of the squeezing pa- 
rameter ^3 versus the parameter 6 for the excitation number 
n = (blue dotted line), n = 10 (orange dashed line) and 
n = 20 (red solid line). The total number of the atoms N are 
chosen as 20. The greatest spin-squeezing is obtained for the 
symmetric Dicke state when the half of the quasi-spins are 
excited (n = iV/2) and 6* = 7r/2. 



FIG. 5: (Color onUne) The polar plot of the squeezing pa- 
rameter ^3 versus the vector difference K. The total number 
of the atoms N ^ the excitation number n and the parameter 
are chosen as 20,4 and 7r/2. Here, ^3 = 1 is the base- 
line to measure whether the dark state is squeezed or not. 
Obviously, the dark state is squeezed only when the vector 
difference K G [-Kc,Kc]. 




FIG. 4: (Color online) The contour plot of the squeezing pa- 
rameter ^3 versus the excitation number n and parameter 0. 
The total number of the atoms N are chosen as 20. The 
greatest spin-squeezing (C3 = 1) is obtained for the symmet- 
ric Dicke state when the half of the quasi-spins are excited 
(n = N/2) and = 7r/2. 

that the squeezing parameter C,\ varies with the excita- 
tion number n for different 9. The dark state is always 
squeezed except without excitation (n = 0) or fully ex- 
cited {n = N). It is noticed that when 7^ 0,7r/2, the 
dark state is squeezed even when the excitation number 
n is equal to the total number of the atoms N. Because 
the dark state actually mix the quasi-spins wave state 



and the photon Fock state, this non-vanishing squeez- 
ing mainly results from the sub-Poisson distribution of 
the photons, which will be explicitly demonstrated in 
the last subsection of the Sec. [Ill Fig. [3] is the polar 
plot of the squeezing parameter p = £,1(0) versus for 
different excitation numbers n. For the state without 
excitation (n = 0), the atomic ensemble always has no 
spin-squeezing because all of the quasi-spins are in their 
ground states. The greatest spin-squeezing is obtained 
for the symmetric Dicke state when the half of the quasi- 
spins are excited {n = N/2) and = 7r/2. Fig. [H is 
the contour plot of the squeezing parameter versus 
the excitation number n and parameter 0. Obviously, 
the greatest spin-squeezing (("3 = 1 drawn by red color) 
can be obtained by adiabatically manipulating the pa- 
rameter ^ to 7r/2 at the half excitation of the quasi-spins 
{n = N/2). 

If the wave vectors between the quantum and the clas- 
sical field are different, the nonzero K results in the non- 
symmetrical Dicke bases in Eq. ([H]) and the spin squeez- 
ing decreases rapidly as the K deviates from 0. Fig. [5] is 
the polar plot of the squeezing parameter p = ^3 (0) ver- 
sus K. The parameters are chosen as N = 20, n = 4, 
and 6 = 7r/2. Since ^3 = 1 is the baseline to measure the 
atomic squeezing, the dark state is squeezed only when 
the vector difference K G [— Kc,Kc], where Kc is shown 
in Fig. [5l Fig. [6] is the contour plot of the squeezing pa- 
rameter (1 versus the vector difference K and parameter 
0. The total number of the atoms N and the excitation 
number n are chosen as 20, 4. In Fig.[6l [— Kc, Kc] is still 
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FIG. 6: (Color online) The contour plot of the squeezing pa- 
rameter ^3 versus the vector difference K and the parameter 0. 
The total number of the atoms N and the excitation number 
n are chosen as 20, 4. Obviously, the dark state is squeezed 
only when the vector difference K G [— Kc, Kc]. 



the boundary for the vector difference K when the dark 
states are squeezed. 



III. INTERNAL ENTANGLEMENT IN DARK 
STATES AND PHOTON STATISTICS 

A. Concurrence in dark states without decoherence 

As the physical observable to quantify the pairwise en- 
tanglement of spin- 1/2, the concurrence is closely related 
to spin squeezing [21]. The concurrence is defined as [2^] 



C — max {0, Ai — A2 — A3 — A4} , 



(32) 



where the quantities Xiii = 1, 2, 3, 4) are the square roots 
of the eigen-values in descending order of the matrix 
product 

Q12 = P12 {(^ly ^ (^2y) PI2 {^ly ^ ^2y) • (33) 

Here, is the two-spin reduced density matrix and 
is the complex conjugate of the pi2' 

For the dark states, the two-quasi-spin reduced density 
matrix f23l 



Pl2 



v+ ii* 

{) w y {) 

y w 

u V- 



(34) 



can be explicitly obtained by tracing out the degrees of 
freedom of all the other quasi-spins. It is clear to verify 
that for the dark states only the elements v±,u,w and y 



of the survive as (0 = K • 1) 

W = w\n_2L + SW, 

a— 2 

u = —wi sin 0, 
y = w cos 0, 



(35a) 

(35b) 

(35c) 
(35d) 



where [23 



n(n — 1) 

N{N -ly 

\{l-2{.l)^{.lal)) 

(n-7V)(n-7V + l) 
N{N -I) ' 

(1 - 



1 

4 
nN 



(36a) 



(36b) 



(36c) 



N{N -I) 

are the corresponding elements of the reduced density 
matrix over the Dicke state \n — N/2) and 

5v+ = B (TV, 0) [n + (TV csc^ 6> + n - 1) A(n, TV, 6)] , 

(37a) 

5v_ = B (TV, 6) [n + (TV cot^ i9 - TV + 2n + 1) A(n, N, 6)] 

(37b) 

6w = B (TV, 0) [-n - {n ^ N cot^ 0) A(n, TV, 0)] (37c) 
with B {N, 6) = cot^ 0/ (TV - 1) and 

(TV + 1) 



A{n,N,0) = 1 



r (n, N, i 



(38) 



(TV-n + 1) 

where T (n, TV, 0) is defined in Eq. (j3Q)) . 

Thus the concurrence of the dark state is given by 

C = max {0, 2 {w |cos ^| - ^/vyl^) } . (39) 

Fig. [3 is the contour plot of the concurrence C versus the 
vector difference K and the parameter 0. In compari- 
son with the one region [— Kc, Kc] for generating atomic 
squeezing, there are three regions for generating pairwise 
entanglement. This is also shown in Fig. [51 which is the 
sections of the atomic squeezing and the concurrence. 
Obviously, in the vicinity of the = tt, the concurrence 
appears but atomic squeezing does not. In Fig. [7] and [HI 
the total number of the atoms TV and the excitation num- 
ber n are chosen as 20, 4. In the next section, the evolve- 
ment of the concurrence as well as the spin-squeezing 
under various decoherence channels will be considered. 



B. Sub-Poisson distribution of the photons 

We have shown that the atomic squeezing for infor- 
mation storage can be generated through the adiabat- 
ically manipulating the angle 0. In another aspect, as 
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FIG. 7: (Color online) The contour plot of the concurrence 
C versus the vector difference K and the parameter 0. The 
total number of the atoms N and the excitation number n are 
chosen as 20 and 4. Here we find three regions for generat- 
ing concurrence rather than one region for generating atomic 
squeezing. 
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FIG. 8: (Color onUne) The sections of the atomic squeezing ^3 
(blue soUd blue) and concurrence C (red dashed line) versus 
the vector difference K. The total number of the atoms iV, 
the excitation number n and the parameter are chosen as 
20,4 and 7r/2. Obviously, in the vicinity of the K = tt, the 
concurrence appears but atomic squeezing does not. 



the photon dressed collective atom state the dark states 
mix the atomic ensemble state and the photon states. 
In our setup, the conservative quantity is the summa- 
tion of the 2;— component collective operator Jz and the 
number operator n = a^a of the photons. The multi- 
partite quantum correlation in dark states actually is the 
component spin squeezing for the dark states, which 
would result from the quantum correlation of the pho- 
tons. Usually, such photonic quantum correlation can be 
measured by the sub-Poisson distribution of the photons 




FIG. 9: (Color online) The sections of the atomic squeezing 
(blue soUd blue) and sub-Poisson distribution Sp (red dashed 
line) versus the parameter 0. The total number of the atoms 
A/", the excitation number n and the vector difference K are 
chosen as 20, 4 and 0. When the increases from to 7r/2, the 
atomic squeezing increases and the sub-Poisson distribution 
decreases simultaneously and vise vesa when the increases 
from 7r/2 to tt. 



defined as 

= 1 ((atataa) - {a^af) . (40) 

If the photonic number fluctuation 

(An2) = (n2)-(n)' (41) 

is smaller than the average of the particle number (n) of 
the photons (sp < 0), the dark states actually has the 
sub-Poisson distribution. 

The straightforward calculation gives 



_ (n-l)n7V^ cot^ i9 

~ (7V-n + 2)(7V-n + l)^ 

iFi (2 - n; A/" - n + 3; -N cot^ 0) 

iFi(-n;7V-n + l;-7Vcot26>) ' 
_ nN cot^ 
~ (7V-n + l)^ 
iFi(l -n;A/'-n + 2;-iVcot^ 0) 

iFi(-n;A/'-n + l;-7Vcot2(9) 



(42a) 



(42b) 



Substituting the relevant expectation values into Eqs. 
(|4Q|) leads to the explicit expression of the sub-Poisson 
distribution of the photons Sp, which always follow the 
atomic squeezing ^3- Fig. [9] illustrates the numerical rela- 
tionship between sub-Poisson distribution and the atomic 
squeezing. When the increases from to 7r/2, the 
atomic squeezing increases and the sub-Poisson distri- 
bution decreases simultaneously and vise vesa when the 
increases from 7r/2 to tt. From the above discussion, 
the greatest atomic squeezing and the concurrence can 
be generated when the parameter 9 is manipulated to 
7r/2, which means all the sub-Poisson distribution of the 
photons converts to the atomic squeezing. 
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In this section, we just consider the adiabatic manip- 
ulation without any decoherence process. In the next 
section, we wih take the typical three atomic decoher- 
ence channels into consideration to estimate the sudden 
death of the atomic squeezing and concurrence. Since 
the coherence time of photons is much larger than atomic 
one, the photonic decoherence such as the dephasing due 
to the imperfect single-photon indistinguishability [24] is 
not considered below. 



From the above expression of the spin-squeezing and 
the concurrence, we notice that if we know the expec- 
tations {Jz) ,(J^)and (crj), and the correlations (ajcr^) , 
all the spin-squeezing and the concurrence can be deter- 
mined. We will give explicit analytical expressions for 
them under three decoherence channels. 



B. Amplitude-damping channel 



IV. ATOMIC SQUEEZING AND PAIRWISE 
ENTANGLEMENT UNDER DECOHERENCE 
CHANNELS 

A. Decoherence channels 

Through the adiabatically manipulating the angle 6>, 
the spin-squeezing and the concurrence are generated for 
information storage. However, due to the existence of the 
environment, they are irreversible depressed and vanish 
eventually. Such decoherence processes are usually de- 
scribed by three typical decoherence channels: the ampli- 
tude damping channel (ADC), the phase damping chan- 
nel (PDC) and the depolarizing channel (DPC) j2i|,[25|. 
They are defined by some maps of the elements of the 
density matrix \i) Since the ADC is an energy- losing 
process and thus it forces the initial state to dissipate 
into the ground state, the map is defined as 

|i)(i|^(l-p)|i)(i|+p|0)(0|, (43a) 

N){jH(l-p)^K)(il,(i7^i). (43b) 

The PDC is phase-losing process and the corresponding 
map is described by 



\i) (j| ^ {l-p + p5i,,)\i) 



(44) 



The DPC is polarization-losing process and the corre- 
sponding map is represented by 



ii) (j|^(i-p)|i)(ii 



(45) 



The behaviors of the spin-squeezing and the concur- 
rence under these decoherence channels will be present 
in the following sub sections. Here, p is the decoherence 
strength whose range is [0,1] . When p = there is no 
decoherence and when p = 1 the decoherence processes 
are completed. 

If we fix the initial values of the spin-squeezing and the 
concurrence of the dark state, the dark state degrades to 
a product state with single quasi-spin state as a compo- 
nent when the decoherence strength p increases from 0. 
Thus the spin-squeezing and the concurrence would have 
sudden death. Before the sudden death happens, the 
information storage can be in progress safely. In other 
word, the longer time the spin-squeezing and the concur- 
rence survive, the better the information storage can be 
achieved. 



Based on the map in Eq. (43) for ADC, one can find 
the following relations for single quasi-spin operators 

{a\{p)) = s{al)^-p, (46a) 
WWAp)) = s^{alal)^-2sp{al)^+p', (46b) 
WWAp)) = s{aial)^{a,/3 = x,y). (46c) 

Hereafter, we define s = 1 — p. Then the corresponding 
relations for the collective operators are 



{J'. (P)> = 

(J' ip)) 



^Jl)^ + {N-l)sp{J,)^ 



s {Jz)o - yP, 

s{3\-sp{j!}^ + iN-l)sp{J,), 
N 



(47a) 
(47b) 



-p{Np + 2s). 



(47c) 



Substituting the relevant expectation values and the cor- 
relation function into Eqs. (|25|26p leads to the explicit 
expression of the spin-squeezing parameters 



es {pY 



^Hp) 



(48) 



with 



(P) = + 2^(^ - {J.)o+Pi^ + s), (49) 

where (^q is the initial spin-squeezing parameter defined 
inEq. ([25]). 

In the same case, in order to investigate the evolution 
of the concurrence, the two-qubit density matrix is shown 
as 



with = 


















yA 








yA 





















= su 


and 



(50) 



A 2 

w = s w 



1 / 1\ 1 

i;sp{a,)^ + -sp, 



(51a) 
(51b) 
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c 



0.01 0.02 0.03 

(b) 




0.64 



,0.05 



0.15 0.3 
P 

(b) 



FIG. 11: (Color online) The sudden death of (a) the atomic 
squeezing C,\ and (b) the concurrence C under PDC. The pa- 
rameters are chosen as N — 20, n — 16 and K — 0. The 
atomic squeezing (^^ always disappear later than the concur- 
rence C. 



FIG. 10: (Color online) The sudden death of the atomic 
squeezing ^3 (blue solid line) and the concurrence C (red 
dashed line) under ADC. (a) Usually, the atomic squeezing 
C3 disappear later than the concurrence C. (b)The atomic 
squeezing ^3 can disappear earlier than the concurrence C. 
(c)The atomic squeezing ("3 may reach its maximum and then 
disappear near the maximum decoherence strength p = 1. 
The parameters are chosen as N — 20, n = 16, = and 
(a)6' = 7r/3,(b)6' = 7r/2,(c)6' = 0.6737r. 



where v±,y^u and w are initial elements defined in Eq. 
(35). Thus we obtain the time evolution of the concur- 
rence as 



C"^ =max{0,C^C2^}, 
= 2 ^s^|cos0| 
= 2 (s^lsin^l - 



w 



(52a) 
(52b) 
(52c) 



Hereafter, we only consider the case at the vector differ- 
ence K = because the quantum correlations decrease 
rapidly when the vector difference deviates from shown 
in Fig. [51 

The subindex A in the above Eqs. (1481501 51.52) rep- 
resents that the evolutions of the spin-squeezing and the 
concurrence are taken under ADC. The sudden death of 
the spin-squeezing and the concurrence implied in Eq. 
(|48l 52) are shown in Fig.[TOl Usually, the atomic squeez- 
ing (1 disappear later than the concurrence C as shown 
in Fig. [TOlfa). However for ADC, the atomic squeez- 
ing (I can disappear earlier than the concurrence C 
(Fig. [TOlf b)). With appropriate parameters, the atomic 
squeezing Q may reach its maximum and then disap- 
pear near the maximum decoherence strength p = 1 
(Fig. [TOlf c)). The parameters are chosen as TV = 20, n = 
16, K = and (a)(9 = 7r/3,(b)l9 = 7r/2,(c)l9 = 0.6737r. 



C. Phase-damping channel 

For PDC described by the map in Eq. ([33]), one can 



find the following relations for single quasi-spin opera- 
tors \JA] 

{^z{p)) = {(Tz)o^ (53a) 

(alalip)) = {alal)^, (53b) 

{aiaj (p)) = {ala})^ (a, f3 = x,y). (53c) 

Then the corresponding relations for the collective oper- 
ators are ( (p)) = {J^}^ , (J^ (p)) = ( J^)^ and 



(J^(P)) 



N 



(54) 



Substituting the relevant expectation values and the cor- 
relation function into Eqs. (|25|26p leads to the explicit 
expression of the spin-squeezing parameters 



4 



(55) 



N' 



In the same sense to investigate the evolution of the con- 
currence, the two-qubit density matrix is present as 



p 
P12 



with relations 







/ 














p 








vl _ 







(56) 



s'^y, and 



s^u. Thus we obtain the evolution of the concurrence as 
= max {0, 2 {s^w |cos 0| - ./vyvZ) ] . (57) 



The subindex P in the above Eqs. (|55|56|57p represents 
that the time evolutions of the spin-squeezing and the 
concurrence are taken under PDC. The sudden death of 
the spin-squeezing and the concurrence implied in Eq. 
(|55|57p are shown in Fig. 11. In contrast of the ADC 
case, the atomic squeezing C,\ always disappear later than 
the concurrence C. 
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0.64 
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0.2 0.4 
P 

(a) 



FIG. 12: (Color online) The sudden death of (a) the atomic 
squeezing (I and (b) the concurrence C under DPC. The pa- 
rameters are chosen as N — 20, n — 4 and K = 0. The atomic 
squeezing (^^ always disappear later than the concurrence C. 



D. Depolarizing channel 

Based on the map in Eq. (|^5]) for DPC, one can find 
the following relations for single quasi-spin operators [14] 



(cTzip)) = s{(Jz)q, (58a) 
{cj14{p)) = s'{aial)^{a,p = x,y,z). (58b) 

Then the corresponding relations for the collective oper- 
ators are {Jz {p)) = s {Jz)q and 



{j'Ap)) = ( J|)o + (1 - s^) 



(J^(P)) 



N 

T' 

3N 



(59a) 
(59b) 



Substituting the relevant expectation values and the cor- 
relation function into Eqs. (|25|26p leads to the explicit 
expression of the spin-squeezing parameters 



N' 



N 



(60) 



In the same sense to investigate the evolution of the con- 
currence, the two-qubit density matrix is present as 



P?2 



with relations 

































_ 






and 



(61) 



5^ + S 



s'^wo 



-v± 

1 



- s 



(62a) 
(62b) 



Thus we obtain the evolution of the concurrence as 

= max jo, 2 (^s^w |cos 0| - ^Jv^v^^ | • (63) 



The subindex D in the above Eqs. (I6UI611 63I55D rep- 
resents that the evolutions of the spin-squeezing and the 
concurrence are taken under DPC. The sudden death of 
the spin-squeezing and the concurrence implied in Eq. 
(|60I63D are shown in Fig. [H In contrast of the ADC 
case, the atomic squeezing (^3 always disappear later than 
the concurrence C. 



V. OPTIMAL TIME FOR GENERATING 
ATOMIC SQUEEZING 

Now, two competitive processes dominate generating 
the spin-squeezing. One is adiabatically manipulating 
the parameter from to 7r/2, which increases the spin- 
squeezing and the concurrence for storing the information 
of phonons into the atomic ensemble. However, the other 
process resulting from the various decoherence channels 
decreases the spin-squeezing and the concurrence, which 
means that the system loses information continuously. 
Therefore, the competition between this two processes 
leads to the existence of an optimal time for storing infor- 
mation, after which the information stored in the atomic 
ensemble is always losing. We can define two time scales: 
if only considering the adiabatic manipulation, one time 
scale is ti representing the time when the spin squeezing 
increases to the half maximum value; if only considering 
the decoherence processes, the other time scale is ^2 rep- 
resenting the time when the spin squeezing decreases to 
the half maximum value. Only when ti > t2 the optimal 
time exists. 

We would like to demonstrate such competition in 
a typical quantum memory based on the cold ^^Rb 
atomic ensemble, which is released from a magneto- 
optical trap at a temperature of about lOOfiK 0]. 
In such quantum memory, the ground state |^), the 
metastable state |m) and the excited state |e) are chosen 
as \^Sii2-,F = l,mi? = l), \^Sii2^F = 2, mi? = — l) and 
|55'i/2, ^ = 2, mi? = 0), respectively. An off-resonant 
cr~ -polarized write pulse contributes to the transition 
from 1^) to |e) and the Stokes photon with a~ polar- 
ization associates with the transition from |e) to |m). 
The dark state for storing the photonic information is 
described in Eq. ([7]). 

To realize the adiabatic creation of the non-classical 
correlation, we adopt the hyperbolic tangent type pulse 
sequence for both Rabi frequencies [26] as 



g{t) = 



1 — tanh( 



t t-T 



1 + tanh(- 



t t-T 



(64a) 
(64b) 



instead of the usual gaussian type ones [27] or solitonary 
type ones [28^. Here, Vt m is the maximum Rabi frequency, 
r is the pulse length and a corresponds to the half width 
of the pulse. The pulse sequences of both Rabi frequen- 
cies are depicted in Fig. [131 Only at the time interval 
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FIG. 13: (Color online) Optical pulse sequences for both Rabi 
frequencies. The red dashed line and black solid line repre- 
sent the Rabi frequencies Q(t) and g{t), respectively. The 
adiabatic following condition is QmT ^ 1. 



[0, r] are both sequences applied to adiabatically manip- 
ulate the dark state in Eq. ([7j). Two Rabi frequencies 
respectively reach their maximum value Qrn at the initial 
time t = and the terminal time t = which guarantees 
that the angle varies from to 7r/2 and thus the infor- 
mation of the photons is stored into the atomic ensemble 
at the time interval [0,r]. 

Actually, due to the adiabatic manipulation, the time 
scale ti for generating spin squeezing is limited by the 
adiabatic condition 



n{t)\dn{t) 



En — Ed 



< 1 



(65) 



with eigenstates \n{t)) (the corresponding eigen-energy 
En) and the time-dependent dark states \dn{t)) (the cor- 
responding eigen-energy E^). Here, 



n{t)\dn{t) 



{n{t)\§-^H{t)\dn{t)) 

En — Ed 



(66) 



describes the variation due to the time-dependent Hamil- 
tonian. Thus the condition for the adiabatic following is 
O^niT ^ 1. The typical pulse length r is about 150//5 [2y|, 
which means the maximum value of the Rabi frequencies 
Vtni should much larger than 6.67 x l{)^Hz. According to 
the pulse sequence in Eq. (|64|) . the time scale ti can be 
determined when (ti) is tuned to 7r/6 as 



2a- 



2r 



(67) 



where r = tanh~^[(6 — 7r)/(6 + 7r)].When a tends to in- 
finity, the longest time ~ l/2r = 75fis is obtained. 

In the other hand, we also assume the exponential de- 
cays for the decoherence processes with the decoherence 
strength 



p{t) = i- 

which leads to the time scale 



t2 = 



-It 



(68) 



(69) 




^C^S 0.2 



FIG. 14: (Color online) The evolution of the atomic squeez- 
ing (blue solid line) and the concurrence C (red dashed 
line). The parameters are chosen as N = 20, n = 4, = 0. 
The maximum of the atomic squeezing and the concurrence 
are respectively obtained at tg ~ 0.127~^ = 120/is and 
tc ^ 0.097"^ 90/is. 



For example, for the cold trapped ^^Rb atomic ensemble, 
the typical dephasing damping rate 7 = 10^ Hz and the 
corresponding dephasing time ^2 = 1ms (tI, f29|. 

The evolution of the spin-squeezing and the concur- 
rence is shown in Fig. [TU Obviously, The maximum of 
the atomic squeezing and the concurrence are respec- 
tively obtained at ^ 0.127~^ = 120 fis and tc ~ 
0.097~^ = 90/i5. Though the Fig. [14] is drawn under 
PDC, the optimal times can also be obtained under ADC 
and DPC. 



Such competition can also be observed by measuring 
the cross-correlation (7|, [S^l 



9s,AS = 1 + CT{t) 



(70) 



as a function of the time, where C is a fitting parameter 
determined by the excitation probability together with 
the background noise in the anti-Stokes channel. Here, 



rW=Kfmax|*(t))r 



(71) 



is the retrieval efficiency, which actually characterizes 
the complete time-dependence properties of the cross- 
correlation. For the dark state we discuss above, to esti- 
mate the adiabatical creation and the decoherence chan- 
nels at the same time, we choose |^max) as the collective 
atomic state \Mn) and thus the retrieval efficiency is ob- 
tained as 



m = 



n\{N-ny. 



2Fi(n- A^,-n;l;e-2T*) 



iFi (-n; N-n + l;-N cot^ e{t)) ' 

(72) 

Here, 2^1 (^^^ 6; c ; z) is generalized Kummer hypergeomet- 
ric function [16|. The retrieval efficiency is depicted in 
Fig. [151 ill which the maximum of the retrieval efficiency 
exists when the same condition ti > t2 is satisfied. 
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ret) 




FIG. 15: (Color online). The retrieval efficiency as a function 
of the time. The parameters are chosen as N — 20, n — 
4, K = 0. The maximum retrieval efficiency is also obtained 
around the 0.127~^ = 120/is. 



Poisson distribution of the photons. The greatest atomic 
squeezing and the concurrence can be generated when 
the parameter is manipulated to 7r/2, which means all 
the sub-Poisson distribution of the photons converts to 
the atomic squeezing. 

We also consider the sudden death for the atomic 
squeezing and the pairwise entanglement under various 
decoherence channels. For the three typical decoherence 
channels, the sudden deaths of the atomic squeezing hap- 
pens later than the sudden death of the concurrence. 

According to the above investigation, an optimal time 
for generating the greatest atomic squeezing and pair- 
wise entanglement is obtained by studying in details the 
competition between the adiabatic creation of quantum 
correlation in the atomic ensemble and the decoherences 
we describe with decoherence channels. 



VI. CONCLUSION 

We investigate the multipartite correlations character- 
ized by the atomic squeezing beyond the pairwise en- 
tanglement in the dark states, which are proposed for 
the quantum information storage based on electromag- 
netically induced transparency mechanism. The atomic 
squeezing and the pairwise entanglement can be created 
by adiabatically manipulating the Rabi frequency of the 
the classical light field on the atomic ensemble. The 
atomic squeezing in dark states converts from the sub- 
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